*********************************************************************************
********		VOU Project - Arrival times & Income 2013 5-year ACS     ********
********		Feb 23, 2021 - A. Waxman								 ********
*********************************************************************************

use ..\..\data\ACS_arrival_time\ACS_Arrival_Time2013,clear

keep if met2013==31080
keep if countyfip==37
keep if empstat==1
drop if incwage<20000
keep if tranwork==10
drop if riders>6
keep if arrives>=300 & arrives<=540
g hrwg=incwage/(2040)
g child=nchild>0
g cpl=carpool>1 | riders>1
twoway (lpolyci hrwg arrives, fc(none) alp(dot) alc(black) clc(black) clw(thick)) ///
	 (lpolyci hrwg arrives if child==1, fc(none) alp(dot) clp(dash)  alc(gs3) clc(gs3) clw(thick)) ///
	(lpolyci hrwg arrives if cpl==1, fc(none) alp(dot) alc(gs8) clc(gs8) clw(thick)) , xtitle("Arrival Time") 	///
	xlabel(300 "5 AM" 360 "6 AM" 420 "7 AM" 480 "8 AM" 540 "9 AM") ///
	ytitle("Hourly Wage in Dollars") plotregion(style(none)) ///
	legend(label(2 "All") label(4 "Children") label(6 "Carpool") ///
	order(2 4 6) row(1) region(lwidth(none)))

graph export ".\results\appendix\figs\Figure5pd.wmf", replace
